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Abstract 



Using a functional-integral approach, we have determined the tempera- 
ture below which cavitation in liquid helium is driven by thermally assisted 
quantum tunneling. For both helium isotopes, we have obtained the crossover 
temperature in the whole range of allowed negative pressures. Our results are 

compatible with recent experimental results on ^He. 
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The possibility of having observed quantum cavitation in superfluid ^He has been first put 
forward by Bahbar and coworkers 0]. These authors have used a hemispherical transducer 
that focusses a sound wave in a small region of a cell where cavitation is induced in liquid 
^He at low temperature. The analysis of their experimental data is complicated by the fact 
that neither the pressure (P) nor the temperature (T) at the focus can be directly measured. 
This makes the determination of the thermal-to-quantum cavitation crossover temperature 
T* to depend on the theoretical equation of state (EOS) near the spinodal point. Using the 
results of Ref. they conclude that T* ~ 200 mK, in agreement with the prediction of 

0. However, using for instance the EOS of Ref. 0], which reproduces the spinodal point 
microscopically calculated by Boronat et al 0,0, the "experimental" result becomes 120 
mK. 

The first detailed description of the cavitation process in liquid helium was provided by 
Lifshitz and Kagan 0, who used the classical capillarity model near the saturation line, and 
a density functional-like description near the spinodal line. More recently, the method has 
been further elaborated by Xiong and Maris [[7|. These authors conclude that there is no 
clear way to interpolate between these two regimes, which makes quite uncertain the range 
of pressures in which each of them is valid. 

In this work, we determine T* for ^He and ^He using a functional-integral approach 
(FIA) in conjunction with a density functional description of liquid helium. The method 
overcomes the conceptual limitations of previous works based on the application of zero- 
temperature multidimensional WKB methods 0, and the technical ones inherent to the use 
of parametrized bubble density profiles 0, thus putting on firmer grounds the theoretical 
results. Moreover, it gives T* in the whole pressure range. 

Thermally assisted quantum tunneling is nowadays well understood (see for example Ref. 
0] and Refs. therein). Let us simply recall that at high temperatures, the cavitation rate, 

1. e., the number of bubbles formed per unit time and volume, is given by 

= Joy e"^^™-/^ , (1) 
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where Aflmax is the barrier height for thermal activation and Jqt is a prefactor which depends 
on the dynamics of the cavitation process. At low T, it becomes 

jQ = JoQe~'-", (2) 

where Smin is the minimum of the imaginary-time action 

S{T) = j^dr JdfC, (3) 

£ being the imaginary-time classical Lagrangian density of the system and the time- 
integration is extended over a period in the potential well obtained by inverting the potential 
barrier. These equations hold provided the rate can be calculated in the semiclassical limit, 
i.e., Smin >> 1; whlch is the present case. For a given value of T, one has to obtain periodic 
solutions to the variational problem embodied in Eq. (j^). Among these many periodic 
solutions, called thermons in Ref. 0, those relevant for the problem of finding T* are the 
ones corresponding to small oscillations around the minimum of the potential, which has an 
energy equal to —AQ^nax- If is the angular frequency of this oscillation, T* = hujp/2Tr. 
It is worth realizing that contrarily to WKB, this procedure permits to go continously from 
one regime to the other: at T*, Eqs. (||) and @) coincide, whereas the WKB approach forces 
to equal a zero-temperature barrier penetrability to a finite-temperature Arrhenius factor 
0J^. Whether this is justified or not, can only be ascertained a posteriori comparing the 
WKB with FIA results. 

To obtain the Lagrangian density C we have resorted to a zero-temperature density 
functional description of the system P,p!0|. This is justified in view of the low-T that are 
expected to come into play (< 200 mK). The critical cavity density profile po{r) is obtained 
solving the Euler-Lagrange equation ||7|,pA| 

where uj{p) is the grand potential density and p is the particle density. AVL^ax is given by 

AVt^ax = J df[uj{po) - Uj{pm)] , (5) 



where pm is the density of the metastable homogeneous hquid. It is now simple to describe 
the dynamics of the cavitation process in the inverted barrier well, whose equilibrium con- 
figuration corresponds to p^ir) and has an energy —/WLmax- We suppose that the collective 
velocity of the fluid associated with the bubble growth is irrotational. This is not a se- 
vere restriction since one expects only radial displacements (spherically symmetric bubbles). 
Introducing the velocity potential field s(r, t), we have 

C = mps — Hi^p, s) , (6) 

where 7i(p, s) is the imaginary-time hamiltonian density. Defining u{r,t) = Vs(r, t), 

n = ]^mp%f - [u}{p) - uj{pm)] ■ (7) 

Hamilton's equations yield 

mp = = —mV{pu) (8) 

OS 

ms = — — . (9) 
dp 

Eq. (g) is the continuity equation. Taking the gradient of Eq. we get 

m— = -V l-mu^ - — } . (10) 
dt [2 6pj ^ ^ 

Thermons p{r,t) are periodic solutions of Eqs. (|^) and (p^O]). From Eq. (^) and using 
Eqs. (^ and (||) we can write 

SminiT) = (j^ dr J df l^mpw^ + Lj{p) - u;(pm)| . (11) 

Within this model, to exactly obtain T* only a linearized version of Eqs. (^ and (ID) around 
Poir) is needed. Defining the T*-thermon as 



p{r,t)=po{r)+p,{r)e''^^\ (12) 

where pi(r) is much smaller than po{r), and keeping only first order terms in u{r,t) and 
pi(r), we get: 




• pi(r) 



) 



(13) 



Here, • pi{r) means that 5uj/5p has to be hnearized, keeping only terms in pi(r) and its 
derivatives. 



shows that its physical solutions have to fulfill p'i(O) = Pi'(O) = 0, and fall exponentially to 
zero at large distances. The linearized continuity equation pi(r) oc —V{pou) imposes the 
integral of pi(r) to yield zero when taken over the whole space. 

We have solved Eq. ([T3|) using seven point Lagrange formulae to discretize the r- 
derivatives together with a standard diagonalization subroutine. The sensibility of the 
solution to the precise value of the r-step has been carefully checked, and in most cases 
a value Ar = 0.25 A has been used. 

For all pressures, only one positive mUp eigenvalue has been found. Fig. 1 (a) and (b) 
shows T* (mK) as a function of P(bar) for "^He and ^He, respectively. In the case of "^He, 
the maximum T* is 238 mK at -8.58 bar, and for ^He it is 146 mK at -2.91 bar. It is worth 
noting that T* is strongly dependent on P in the spinodal region, falling to zero at the 
spinodal point (see also Ref. 

We display in Fig. 2 the pi(r)-component of the thermon ([121) in the case of ^He (a 
similar figure could be drawn for ^He). For large bubbles, pi(r) is localized at the surface: 
the thermon is a well defined surface excitation. It justifies the use of the capillarity approx- 
imation near saturation, or more elaborated approaches, like that of Ref. [§], that consists 
in a simplified one-dimensional model in which the oscillations are just described by rigid 
displacements of the critical bubble surface. 

When the density inside the bubble becomes sizeable, a mixed surface-volume thermon 
develops, which eventually becomes a pure volume mode in the spinodal region. This mode 
can no longer be described as a rigid density displacement, and the above mentioned models 
fail: the exact T* is higher than the prediction of the rigid surface displacement model 
because volume modes involve higher frequencies. 



Eq. (|I^) is a fourth-order linear differential, eigenvalue equation. A careful analysis 
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To determine which of the T*(P) shown in Fig. 1 corresponds to the actual experimental 
conditions, we have calculated the homogeneous cavitation pressure 0,0] • It is the one 
the system can sustain before bubbles nucleate at an appreciable rate. We have solved the 
equation 

1 = [Vt), J (14) 

taking J=Jt and 

Vq = AnRl/S represents the volume of the critical bubble, for which we have taken Rc =10 
A. For T < T*, Jt has to be replaced by Jq. Lacking of a better choice, we have taken 
Jqq = Jor(T=T*), and for the experimental factor (Vt)e (experimental volumextime), two 
values at the limits of the experimental range [|l],0], namely 10^^ and 10^ s. For ^He it 
yields P/i=-8.57 bar and -8.99 bar, respectively. The corresponding values for ^He are -2.97 
and -3.06 bar. This means that for both isotopes Ph is close to the spinodal pressure. Table 
1 displays the associated T*-values. 

The crossover temperatures are similar to those given in Ref. 0, although different 
functionals have been used in both calculations. As a matter of fact, this is irrelevant, 
since both functionals reproduce equally well the experimental quantities pertinent to the 
description of the cavitation process. 

An explanation for the agreement between these calculations can be found in Ref. 0. In 
that work, using a simplified one-dimensional model in which the oscillations were modelled 
by rigid displacements of the bubble surface, the cavitation process was described within FIA 
from T=0 to the thermal regime. It was shown that thermally assisted quantum cavitation 
only adds small corrections to the T=0 "instanton" solution (formally equivalent to WKB 
if Smin >> 1) in the quantum-to-thermal transition region. 

Let us recall that the formalism used in Ref. |^ to describe quantum cavitation is a 
multidimensional WKB one, appropiated for a T=0, pure quantum state with a well defined 
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energy value. This approximation is well known to fail for energies close to the top of the 
barrier. On the contrary, the FIA here adopted deals with thermally mixed quantum states, 
making it possible to smoothly connect quantum and thermal regimes Q. Besides, it is 
technically complicated to obtain the E=0 instanton solution to Eqs. (§) and (|T0|) without 
using some numerical approximations that might be unworkable in more complex physical 
situations, like that of a "^He-^He liquid mixture. We also want to stress again that, to 
determine the quantity of experimental significance, namely T*, only the thermon solution 
of the much simpler eigenvalue Eq. (|13|) is required. 

To conclude, within density functional theory, we have performed a thorough description 
of the quantum-to-thermal transition in the process of cavitation in liquid helium based 
on the functional- integral approach. Our quantitative results (see also Ref. 0) indicate 
that the crossover temperature is below 240 mK for '^He, and below 150 mK for ^He. The 
experiments on ^He yield results which, depending on which equation of state is used, are 
in the 120-200 mK range. Given the present uncertainties in theoretical and experimental 
results as well, we consider the agreement as satisfactory. 
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tine) Grant No. PID 97/93 and by the IN2P3-C1CYT agreement. 
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FIGURES 

FIG. 1. (a) T*(mK) as a function of P(bar) for ^He. (b) Same as (a) for ^He. 

FIG. 2. Referring to ^He, we show: (a) the particle density profile po{r) (solid line) and the 
density pi(r) (dashed line) for P=-4.59 bar. (b) Same as (a) for P=-8.35 bar. (c) Same as (a) for 
P=-9.16 bar. pi(r) is drawn in arbitrary units, po{r) in and r in A. 
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TABLES 

TABLE I. Crossover temperatures for two different values of the experimental volume times 
time. 
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